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Abstract 

By  using  a  variational  principle  for  scattering  by  rough  interfaces,  the  functional 
derivative  of  the  scattered  part  of  an  acoustic  field  originating  from  a  point  source 
is  derived  directly  in  the  time  domain.  This  general  result  is  specialised  to  the  case 
of  scattering  by  surfaces  on  which  fields  satisfy  Dirichlet  boundary  conditions.  It  is 
shown  that  when  the  small-  slope  appoximation  is  used  for  surface  fields,  a  simple 
result  can  be  obtained  for  monostatic  scattering. 
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The  small-slope  approximation  was  developed  by  Voronovich  to  treat  scattering  of  plane,  t  Specie 
monochromatic,  scalar  waves  by  rough  interfaces  [1].  Numerical  experiments  have  shown  | 

that  the  zeroth  order  small-slope  approximation,  generally  out-performs  the  conventional  •' 

KirchhofF  approximation  without  added  computational  work  [2].  The  point  of  the  present 
work  is  to  show  how  the  lowest  order  small-slope  approximation  is  expressed  when  point, 
pulsed  sources  are  considered.  It  turns  out  that  in  the  case  of  monostatic  scattering  from  a 
rough  surface  on  which  Dirichlet  boundary  conditions  are  satisfied,  the  zeroth  order  small- 
slope  approximation  takes  on  a  particularly  simple  form.  This  result  is  readily  compared 
to  a  result  for  monostatic  scattering  in  the  KirchhofF  approximation,  Eqs.(31,  32)  below. 

The  latter  was  developed  by  Berry  [3]  for  the  purpose  of  inferring  surface  roughness  from 
the  shape  of  pulse  echoes. 

Recently  Dashen  and  Wurmser  [4]  rediscovered  a  perturbation  scheme  for  acoustic 
and  electromagnetic  scattering  which  was  previously  developed  by  Burrows  [5]  for  elec¬ 
tromagnetic  scattering  by  surfaces.  The  small-slope  and  other  approximation  schemes 
are  conveniently  derived  from  this  perturbation  theory,  at  least  for  plane  wave  scattering 
[4,  6].  Rather  than  synthesize  plane  waves  to  obtain  the  present  results  for  point  sources, 
the  perturbation  result  of  Dashen  and  Burrows  will  be  rederived  directly  for  time-domain 
Green’s  functions.  In  the  plane  wave  case,  once  a  formula  is  obtained  for  the  derivative 
of  the  scattered  field  with  respect  to  upward  rigid  translations  of  the  surface,  it  is  easy 
to  derive  a  result  for  the  scattering  amplitude.  This  is  because  the  plane  wave  scatter¬ 
ing  amplitude  changes  by  a  simple  phase  factor  when  the  scattering  surface  is  rigidly 
translated  [6].  There  is  no  such  simple  way  to  obtain  comparable  formulas  for  Green’s 
functions  since  each  plane  wave  making  up  the  Green’s  function  changes  phase  by  a  differ¬ 
ent  amount.  However,  as  will  be  shown,  for  the  case  of  monostatic  scattering  it  is  possible 
to  integrate  the  functional  derivative  of  the  field  with  respect  to  the  surface  shape  after  . 
the  small-slope  approximation  is  made  for  the  surface  fields.  This  proceedure  results  in  a 
simple  expression  for  the  echo  of  a  pulse  heard  at  the  location  of  the  origin  of  the  pulse. 

In  the  Burrows-Dashen  result,  the  perturbation  of  the  scattered  field  is  expressed  in 
terms  of  products  of  the  total  field  or  its  derivatives  on  the  scattering  surface.  In  contrast, 
the  Helmholtz- KirchhofF  formula  expresses  the  scattered  field  itself  in  terms  of  products 
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of  unscattered  and  total  fields  on  the  scattering  surface.  In  the  Burrows-Dashen  result, 
reciprocity  is  manifest,  but,  in  general,  only  derivatives  of  scattered  fields  are  obtained. 
In  the  Helmholtz- Kirchhoff  result,  reciprocity,  though  respected,  is  not  manifest. 

In  section  II  the  Burrows-Dashen  perturbation  result  will  be  rederived  directly  for 
point  sources  in  the  time  domain  .  In  Section  III  the  new  result  for  monostatic  scattering 
will  be  derived  from  the  perturbation  result.  It  will  be  compared  to  a  corrected  version  of 
Berry’s  result  based  on  the  Kirchhoff  approximation. 

II  The  Variational  Principle  and  Perturbation  Theory 

The  Burrows-Dashen  perturbation  scheme  can  be  derived  from  the  principle  of  least  action. 
Kohn  [7]  and  later  Altshuler  [8]  used  specializations  of  the  least  action  principle  to  derive 
variational  principles  for  scattering  amplitudes  and  the  field  scattered  by  plane  waves  for 
the  case  of  quantum  potential  scattering.  Here  these  results  are  generalized  to  the  case 
of  point  sources  and  point  receivers  in  the  time  domain  in  media  with  variable  densities 
and  sound  speeds.  Furthermore,  layered  media  with  irregular  fluid-fluid  interfaces  (  across 
which  pressure  and  normal  velocity  are  continuous)  will  be  considered.  To  treat  point 
sources  it  will  be  assumed  that  trial  functions  vanish  at  infinity  because  there  is  a  small 
dissipation.  The  general  variational  principle,  which  is  just  the  expression  of  the  principle 
of  least  action,  is  interesting  because  it  embodies  the  applications  of  Green’s  theorem  that 
go  into  a  variety  of  approximations  [9].  However,  it  will  be  used  here  only  to  treat  the 
very  special  case  of  a  single  Dirichlet  boundary,  constant  sound  speed  and  density,  and 
monostatic  geometry.  One  advantage  of  approaching  the  perturbation  formula  through 
the  variational  principle  is  that  some  of  the  complications  of  continuing  normal  derivatives 
away  from  the  unperturbed  surface  which  were  encountered  by  D ashen  and  Wurmser,  can 
be  avoided.  It  should  also  be  pointed  out  that  formulas  similar  to  those  presented  by 
Dashen  and  Wurmser  not  only  were  given  by  Burrows,  but  also  by  Garabedjian  in  1955 
and  apparently  by  Hadamard  even  earlier  [10]. 

II.  1  The  Variational  Principle 

Consider  acoustic  waves  in  a  medium  with  spatially  variable  ambient  density,  p( r)  and 
sound  speed,  c(r).  The  Green’s  function  Ge,  (e  for  exact)  for  the  pressure  at  a  point  r'  at 
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time  t'  arising  from  an  impulsive  source  at  r0  at  time  ta  satisfies  the  wave  equation 

p(r')V'^V'Ge(r'’*';r<”*o)  "  ^^(I'.t'jr.,*.)  =  S(t'  -  r 0)6(t'  -  t0).  (1) 

Boundary  conditions  are  built  into  this  equation:  Ge  and  n  •  VGe/p  must  be  continuous. 
If  they  did  exhibit  jump  discontinuities,  the  spatial  gradient  operators  would  give  rise  to 
6  functions  which  appear  nowhere  else  in  Eq.l.  These  boundary  conditions  express  the 
requirement  that  the  pressure  G  and  the  normal  velocity,  (n  •  9tVG)/p,  be  continuous. 
Causal  Green’s  functions  are  assumed  here,  so  that 


Ge{r  it  ;r0,t0)  —  0 


if  t'  <t0.  Reciprocity  in  the  time  domain  is  expressed  by  [11] 

-^-Ge(r',  t r0 ,  t0)  =  -r-rGe( r0, -t0\  r\ -tf).  (2) 

P(rj  p(r0) 

A  variational  estimator  T  for  Ge  is  a  functional  of  trial  Green’s  functions,  G\  and  G2, 
T(Gi,G2),  such  that  if  G,  =  Ge  +  (gi,  t  =  1,2,  (where  gi  is  arbitrary  except  for  certain 
boundary  conditions)  then  T (Cn,G2)  -  Ge  =  0(e2)  for  small  c  [12].  Consider  two  fluids 
occupying  regions  fl+  and  ft-  which  are  separated  by  an  interface  S.  Within  each  region 
the  sound  speed  and  density  are  continuous,  but  across  the  interface  they  are  allowed  to 
have  jump  discontinuities.  Trial  functions  and  the  exact  Green’s  functions  are  assumed  to 
vanish  at  infinity,  since  a  small  attentuation  is  assumed.  Trial  functions  are  assumed  to  be 
causal  as  well.  In  the  following  equation  T  is  a  variational  estimator  of  the  Green’s  function 
pGei  which  depends  parametrically  on  the  source  location,  (r*,,^)  and  the  receiver  (or  / 
inal)  location  (r /,!/). 


^(Gjt  T/,  t/i  Gi,  t07 10)  —  ^  ^/)  4"  j 1/i  1#) 

+  f"  df  t  dr'IV'Gjfr',  r,,  -(,)  •  -i-V'G.fr',  <■,  r„,  1.) 

J — oo  J O+uO-  P\*  ) 

-dt'G2(r’,  -if) lfr't «bt *o)] 

+  dt'  fs  dS[J(G2(r', rh  •  V'G1(i',«';  t0,t0)) 
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+  A[W  •  V'G2(r', -t';  r /,  J[Gi(r\  t';  r0,  <„)]],  (3) 

where  n  is  the  unit  normal  to  the  interface  5  from  ft-  to  fi+.  The  operator  J  here  is  the 
jump  across  the  interface  S : 


J[/W]  =  /(r+)-/(r-). 

Likewise,  the  averaging  operator  A  is 

4/«]  =  [/(r+)  +  /(r')]/2. 

When  G\  =  G2  =  Gtriah  the  volume  integral  in  Eq.(3)  in  just  the  negative  of  the  Lan- 
grangian  for  free  scalar  fields  obeying  the  Helmholtz  equation  [14].  The  linear  terms  in 
Gi  and  Gi  account  for  the  fields  here  being  driven  by  ^-function  sources,  and  the  surface 
integrals  are  designed  to  account  for  boundary  conditions  at  an  interface 

To  show  that  T  is  a  variational  estimator  of  Gt/p ,  write  Gi  =  Ge  +  egi ,  *  =  1,2.  Then 
integrate  by  parts,  separately  over  ft*  and  over  time,  to  remove  the  gradient  operation 
and  time  derivatives  from  gi  where  possible,  to  show 

r(<7.,r/,*/;<7«,r<>,«0)  =  ,  f/j  r«,  f»), 

^/j fyi gu toi =  ,  \Gt{rj,tf',T0,t0),  (5) 

PKTf) 

and 

r(<?2>  ^e»  *"o»  fo)  =  /  \  ~foi  T/,  — t/).  (6) 

P\To) 

This  means  that 

r(G2,  tj,  tj\ G\,  t0,  t0)  =  ^  jGe(r/,ty;  r0,<()  +  cgi ,  r0,  t0),  (7) 

where  f  is  the  bilinear  form  T  -Gi/p-  Gi/p.  Since  f  is  bilinear,  the  last  term  in  Eq.(7) 
is  order  e2,  and  Eq.(7)  shows  that  T  is  a  variational  estimator.  Note  that  if  the  Gi  are 
replaced  by  a  single  function  ip,  the  variational  estimator  T  is  the  classical  action  for  the 
field  driven  by  ^-function  sources.  If  ip  is  an  arbitrary  snperpostion  of  two  fields,  G\ 
and  <?2,  the  variational  estimator  is  seen  to  follow  from  the  principle  least  action.  The 
integral  along  the  surface  5  in  the  action  represents  a  singular  part  of  the  Lagrange  density. 
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Physically  it  arises  because  there  is  work  done  if  the  interface  moves  under  the  influence 
of  a  pressure  difference  across  the  interface. 

For  the  case  of  a  Dirichlet  surface  replace 

•/[G<]  -  G, 


and 

y4[n  •  VGj/p]  -♦  n  -VGi/p, 

and  for  the  Neumann  case  drop  the  surface  integral  over  S  entirely.  Note  that  in  the 
integrals  G2  and  Gi  always  appear  with  the  same  arguments.  Hence  one  could  omit  the 
arguments  without  danger  of  confusion.  Note  too,  that  any  number  of  regions  ft  could  be 
considered  as  long  as  interface  integrals  are  included  in  the  variational  estimator. 


II.  2  Perturbations 


The  change  in  the  Green’s  function  caused  by  a  small  change  in  the  shape  or  location 
of  the  interface  5  is  found  from  the  variational  estimator  by  choosing  the  trial  functions 
Gi  and  Gj  to  be  the  Green’s  functions  associated  with  a  neighboring  interface,  So-  It  is 
assumed  that  the  densities  and  sound  speeds,  p*(r)  and  c*(r)  can  be  continued  smoothly 
from  ft*  to  ft*  along  with  the  Green’s  function  on  either  side  of  the  interface,  G*  [13]. 
Then  the  variational  estimator  shows  that  to  first  order  in  the  displacement  of  the  surface 
from  So  to  S,  the  Green’s  function  G  associated  with  S  is  given  by 


^  jG(r/,t/;r®,fa)  ss  ^  ^C?a(r0,  t0]Tf,  f/)d*  *<>»$•) 


+  dtf  f  dr'[V'G0(r',  r,,  -tf)  -  V'G0(  r\  *  r0,  t0) 

J- 00  J  O+uO-  P\r) 

-dt'Go(r\  —t'\  r  /, -tf)  dt>G0W,  t'; r0,  f0)] 

+  f  M  Js d5[J[Go(r,,-t/;r/,-t/)]A[~^n,.V,G0(i',<,;r0,ttf)] 

+  A[n'  •  V'G0(r',  r ,,  -tf)-^}J[Go(r',  t';  re,  t0)]].  (8) 
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On  the  other  hand,  the  unperturbed  Green’s  function,  Go,  is  exactly  given  by  the  same 
formula  with  the  regions  ft*  replaced  by  the  regions  ft*.  The  interface  integral  vanishes 
because  Go  is  continuous  across  So,  by  assumption.  If  this  analogous  expression  for  Go  is 
subtracted  from  Eq.8,  the  left  side  is  the  perturbation  of  the  Green’s  function.  The  linear 
terms  cancel  on  the  right.  This  leaves  the  difference  of  the  volume  integrals  which  is  of 
order  of  the  volume  between  the  interfaces  and  the  surface  integral  which  is  first  order 
because  the  jump  in  Go  vanishes  on  the  neighboring  surface  So-  Separating  the  integration 
over  the  upper  +  volumes  from  the  -  volumes  gives 


P(«7) 


6G(rf,tf;r0,t0) » 


+ £  *  L  -  4  *‘G°  -  d'G°wmd'G^ 

+ £  *  L  -  4  ^ 

+  jH  H  J  <fS[./(Go|A[^pjll'  ■  V'G„)]  +  A[n'  ■  V'Go^j]J[Go]]. 


(») 


In  this  equation,  the  arguments  of  the  Green’s  functions  have  been  omitted;  they  axe  to 
be  understood  as  they  appear  in  Eq.8.  If  the  surface  5  is  obtained  from  So  by  normal 
displacement  £(r)  into  ftj,  the  volume  ftp  is  diminished  by  /  dS£  and  the  volume  ftp  is 
increased  by  this  amount.  The  jump  in  the  surface  integral  is  found  from 


Gt  ~  Go  Is  »  G+  -  Go  Isb  +  f(n  •  VG+  -  n  •  VGp  )|*  =  Z*  *  VG+(1  -  fo.  (10) 

Thus,  using  the  boundary  conditions  satisfied  by  Go,  the  perturbation  of  the  Green’s 
function  satisfies, 

^£y^G(r/,t/;r0,t0)» 

-  £ 4  “  w'c°+  •  -  v'°”'  •  mfi  V'G;1 

+ £ d{  L  ds&r'M’Gz^w^a‘<  - 

+  r  dt'  /  dSt( V)2n  •  VG+-£n  •  VG+(1  -  £).  (11) 

J- 00  JSo  p  P 


6 


To  first  order  it  doesn’t  matter  if  the  integrations  are  taken  over  5  or  Sq.  This  result  can 
be  simplified  still  further  by  using 

VGo  •  V<?o  =  n  •  VGon  •  VGo  +  V||Go  •  V||Go 

and  the  boundary  conditions,  which  require  that  V||Go  and  dtGo  be  continuous.  Further¬ 
more,  for  infinite  interfaces,  as  in  layered  media,  it  is  convenient  to  refer  the  displacement 
of  the  surface  to  a  horizontal  surface.  Then  the  displacment  is  described  by  its  projection 
onto  the  vertical  axis,  6h( R).  Surface  areas  transform  by  dS£  =  dRJih.  The  perturbation 
of  the  Green’s  function  is  now  given  by 

+ it'  J  iR'ik(r')  •  VGJ ■  VGJ (1  -  £)  -  VnG„  •  V|,G„(i  - 

+8,.Gj-a,.Gj  (<;1(r,+w^+)  -  '  (12) 

The  analogue  of  this  result  for  plane  wave  scattering  can  be  found  in  Re£s.[4,  6].  In  a 
layered  medium  it  is  clear  that  the  perturbation  of  the  Green’s  function  will  be  a  sum 
of  contributions  from  the  perturbation  of  each  interface  similar  to  that  given  above.  To 
lowest  order  these  perturbations  simply  add  together.  In  the  case  of  Dirichlet  surfaces  only 
the  first  term  of  the  integrand,  involving  normal  derivatives,  contributes  with  p~  =  0. 

It  should  be  noted  that  the  surface  So  is  not  generally  a  flat  surface  and  that  as  a 
consequence  Go  is  not  easily  calculated.  Nevertheless,  Eq.(12)  is  useful  in  that  it  provides 
an  exact,  albeit  formal,  expression  for  the  functional  derivative  of  the  Green’s  function 
with  respect  to  the  surface  shape,  h( r),  since  Eq.(12)  is  exact  to  first  order  in  Sh. 
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Ill  The  Small-Slope  Approximation 

Consider  now  the  special  case  of  a  Dirichlet  surface.  Let  the  scattering  surface  be  described 

by 

z  =  <h(R), 

where,  as  above,  R  is  a  two-dimensional  vector  in  the  x  —  y  plane.  Denote  the  Green’s 
function  associated  with  this  surface  by  Gc  and  consider  what  happens  when  the  surface 
is  perturbed  from  z  —  eh( R)  to  z  —  (e  +  de)h(R).  Then  6h  is  given  by 


6h(  R)  =  {de)h(  R). 


(13) 


Dividing  the  perturbation  of  Gt  by  de  gives  the  following  exact  (because  Eq.(12)  is  exact 
to  first  order  in  Sh )  result  for  the  derivative  of  the  Green’s  function  with  respect  to  e. 

dG 

(r/»  f/i  *o,  fo)  = 

J°°  dt'J dR'^ROn,  •  VG(r/,t/;r',t,|{€/i})iif  •  VG(r',t/;r0,t0|{ch}).  (14) 

On  th^  right  side  of  this  equation,  the  functional  dependence  of  the  Green’s  functions  on 
the  scattering  surface  eh  has  been  made  explicit.  The  Green’s  function  for  the  desired 
surface,  h( R),  namely  G\,  is  found  by  integrating  from  e  =  0  to  e  =  1.  The  function  Gt=o 
is  just  the  flat  surface  ( z  =  0)  Green’s  function  which  can  be  found  fairly  simply  if  the 
density  and  sound  speed  depend  only  on  the  depth  z.  In  performing  the  integration  over 
e,  assume  that  the  order  of  the  integration  can  be  changed  as  needed,  and  integrate  over 
e  first.  Then  change  the  variable  of  integration  from  e  to 

z'  =  eh(  R').  (15) 

One  then  obtains  for  the  Green’s  function  associated  with  the  surface  h( R), 

G(r/,*/;r0,t<,|/i)  -  G(r/,t/;ro,*o|0)  = 

/oo  t  t  h(R!)  j 

dtj  dR'JQ  n(.VG(r/,t/;4<'|{^MK  VG(r',l,;r0,t0|{^M).  (16) 

In  these  expressions,  the  normal  to  the  surface,  n<  is  given  by 


z  -  fVf>(R')  ih(RQ  -  yVhfR') 

“  y/1  +  (cV/»( R'))’  ~  v/h(R'1)2  +  (^(R'))3’ 


(17) 
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In  addition,  in  the  last  equation,  the  intermediate  vector,  r\  is  understood  to  be  given  by 


r'  =  (R',-0- 


(18) 


It  should  be  emphasized  that  Eq.(16)  is  an  exact  result  for  Dirichlet  surfaces  even  when 
there  are  sound  speed  and  density  variations. 

Now  assume  that,  in  fact,  sound  speed  and  density  are  constant.  The  small-slope 
approximation  consists  of  replacing  the  normal  gradient  of  the  Green’s  function,  Gt,  by 
twice  the  vertical  derivative  of  the  free-space  Green’s  function,  G /ree,  the  Green’s  function 
in  the  absence  of  any  boundaries.  In  the  present  case,  G /ree  is  given  by 


G/ree(rf,tj;r0,t0)  =  -6(tf  -t0-  |r/  -  r0|/c)/(4ir|r/  -  rc|).  (19) 


Note  that  this  differs  by  a  sign  and  a  factor  of  4t  from  the  conventional  Green’s  function, 
because,  for  the  purposes  of  deriving  the  variational  principle,  it  was  convenient  to  have 
the  wave  equation  driven  by  a  simple  6-function.  It  is  now  convenient  to  consider  the  field 
V>(r,  t)  produced  by  a  pulse  F(t)  originating  at  r0.  In  the  small-slope  approximation  the 
difference  between  a  pulse  received  when  there  is  a  rough  surface  and  the  pulse  received 
when  there  is  a  flat  surface  located  at  z  =  0  is  given  by 

V>(r /,*/)  -  i>]iat{rj,tj)  =  /  *  dt0[G{T/,tj\T0,t0\h)  -  G(rf,tf;ro,to\0)]F(to) 

«/  — OO 


4 

16jt2 


h'l 


mR') 


Jja  q  -  *o  -  (|r/  -  r'l  +  Ir'  -  r0|)/c 

dZ  ^ - kF^-rol - • 


(20) 


The  point  of  this  section  is  to  show  that  the  z1  integral  can  be  performed  when  the  source 
and  receiver  coincide,  i.e.,  when  r/  =  r0  =  r.  The  main  result  is  given  in  Eq.31  below. 
It  is  a  relatively  simple  expression  but  there  is  a  good  deal  of  algebra  between  here  and 
there.  To  see  how  this  works  out,  first  carry  out  the  differentiation  and  then  set  the  source 
position  equal  to  the  receiver  postion.  Let 


P 


|r- 


Without  loss  of  generality,  suppose  that  r  =  (0,  z).  Then 

1  t  /MR') 

^(r/»</)  “  J  dR'  dJ 
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(22) 


cV2 


cp3  p* 

To  deal  with  this  integral  use  the  identities  [3] 

pc  R'  •  VF 


F;  =  - 


and 


F  =  - 


2  R'2 

pcR'VF 


(23) 


Let  /i  be  the  integral  of  the  term  involving  the  second  derivative 

h'l 


,  (h{R,)dJF"(t-2p/c)(z-z')> 


c'p1 


hi 


/  fh(R) 


o  2pc  R'2  p2 


Write  the  last  integrand  as  a  divergence  in  two  dimensions  and  use  the  identity 


V  •  ^  =  2»*(R). 


(24) 


(25) 


When  the  divergence  operator  is  passed  outside  the  inner  integral,  a  term  in  Vh  is  produced 
from  the  upper  limit  of  the  inner  integral.  In  this  way  I\  can  be  shown  to  be  given  by 


h 


dR 


,R'  •  Vh(R')  F  (2  -  h(R'))2 


M2 


2pc 


2s-  f*<°>  j  /F(t-2(2-z/)/c) 


2*  p  , 

+  4s-2  y0  2  2c(z-2/) 

4t2  y  y0  2cp3  p3 


(26) 


The  last  term  in  this  expression  can  be  combined  with  the  second  integrand  in  Eq.(21)  so 
that  the  scattered  field  at  the  receiver  is  now  given  by 


„  ,  ,  ,  ,1  /^R'-v/KRO  r  (2-h(R#))3 

-  tWr/.t,)  =  _/  JR' - p - ^ ^ - • 

2»  /Wl  ,  ,T(I  -  3(;  -  xQ/t) 

+4  W„  2c(j-S) 


+ 


^SiK 'I 


,  fh{R>)  dz,F(t-2p/c)(z-Z')* 


(27) 


The  next  to  last  term,  which  will  be  denoted  by  I2,  can  be  treated  similarly  to  I\,  with 
the  result 

h  =  ThJd* 


t,R'  •  Vh(R')  F(t-2p/c)(z-h(R'))> 


R' 

2  it  /-M°)  ,F(t-2{z-z')/c) 
+TT-7  dz 
Jo 


16tt2 


(z  -  z')! 

(28) 

4?r2  ^  y0  p4  p2 

It  seems  a  bit  magical,  but  when  this  expression  for  I2  is  used  in  Eq.(27),  the  integrands 
involving  F  cancel  one  another.  Furthermore,  since 

v>  r 

(29) 


d  F(t  -  2(z  -  z')/c)  _  2 F‘  __F_ 

dz'  z  —  zf  c(z  —  z')  +  z  —  z1' 

the  integrals  over  z'  alone  can  be  performed.  Finally,  the  scattered  pulse  is  given  by 


^(r/,t/)-t/»//o«(r/,t/)=  ^ 


F(i  -  2 (z  -  h( 0))/c)  F(t  -  2 z/c) 

2 (z  -  MO)) 


F(t-2p/c)  j  F’(t-2p/c) 


2z 

(*  -  MRQ)2 

p2 


(30) 


4p2  2 cp2 

Bee?  ^e  the  sign  of  the  free-space  Green’s  function  is  negative,  the  pulse  reflected  from  a 
flat  surface  at  z  =  0  is  just 


+ 


1  F(t-2z/c) 


4?r  2  z 

This  means  that  in  the  small-slope  approximation,  the  scattered  portion  of  the  pulse 


0S54  =i>  ~  direct 


is  given  by 


rl>SSA(rf,tf)  = 


1  F(t  -  2(z  -  M0))/c) 


4 -k  2 (z  -  MO)) 

The  only  approximation  used  in  obtaining  this  result  is  the  small-slope  approximation. 

A  similar  result  can  be  obtained  for  the  Kirchhoff  approximation  by  using  the  same 
sorts  of  manipulations  described  here: 

.Kir,  }  = 

*  {/,j)  4*  2(z  —  MO)) 
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+ 


4  T2 


r  R'.vft(R') 

\F(t-2p/c)  F'(t-2p/c)] 

J  dR  R'2 

4/>2  2  cp2 

(32) 


This  last  result  differs  slightly  from  Berry’s  expression  [3]  in  that  no  far-field  assumptions 
have  been  invoked. 

Both  the  Kirchhoff  and  the  small-slope  approximations  give  the  exact  result  for  flat 
surfaces  through  z  =  h(Q).  Both  show  thaf  scattering  takes  place  only  because  the  gradient 
of  the  surface  shape  is  non-vanishing.  The  two  approximations  differ  only  in  that  the 
small-slope  approximation  includes  the  factor 

(*  ~  MR'))2 


This  factor  is  the  sine  squared  of  the  local  angle  of  propagation  at  the  surface.  For  plane 
waves,  the  Kirchhoff-  and  Small-slope  approximations  differ  by  a  factor  of  kzqz/k 2  [6]. 
In  backscattering,  vertical  components  of  incident  and  out-going  wavenumbers  are  equal. 
Thus  the  difference  exhibited  here  between  the  Kirchhoff  and  Small-slope  approximations 
for  point  sources  is  consistent  with  the  difference  observed  for  plane  waves. 


IV  Discussion 

The  primary  results  presented  above  are  the  variational  principle  and  the  subsequent 
derivation  of  what  is,  in  effect,  an  exact  expression. for  the  functional  derivative  of  the 
scattered  field  with  respect  to  the  surface  shape.  This  was  used  in  the  previous  section 
to  derive  a  simple  expression  for  the  small-slope  approximation  for  point  sources  in  the 
time-domain.  The  only  approximation  used  to  derive  Eq.(31)  is  Eq.(19).  The  variational 
principle  should  have  general  utility.  It  was  shown  in  Ref.[9]  that  a  variety  of  approxi¬ 
mations  can  be  based  on  the  variational  principle.  If  the  trial  Green’s  functions  are  the 
same,  the  variational  expression  is  guaranteed  to  be  reciprocal.  Furthermore,  scattering 
is  described  even  if  the  trial  Green’s  functions  don’t  exhibit  scattering.  For  example,  one 
might  try  to  describe  reverberation  in  the  ocean  by  using  trial  Green’s  functions  calculated 
by  the  method  of  adiabatic  normal  modes. 

What  is  surprising  and  annoying  is  that  so  much  algebra  is  required  to  demonstrate  the 
small-slope  approximation.  Such  a  simple  result  deserves  a  simpler  derivation.  Further¬ 
more,  this  result  does  not  seem  to  be  generalizable  to  bistatic  scattering,  nor  to  variable 
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sound  speeds,  at  least  not  without  further  approximations.  It  may  represent  one  of  a  very 
few  situations  in  which  there  is  something  simple  about  the  special  case  of  backscattering. 
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